Magnetic properties of disordered CoCu alloys : A first principles approach 
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Crystalline Coi,Cui_a; alloys show interesting magnetic behavior over the entire concentration 
regime. We here present a fully self-consistent first principles electronic structure studies of the 
electronic structure and magnetic properties of the system. We present results for the variation of 
density of states, magnetic moment, spin susceptibility and Curie temperature. 
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I. INTRODUCTION 

CoCu alloys have been extensively studied experimen- 
tally. Earlier investigations, particularly those of the 
magnetic properties, include high-temperature magnetic 
moment for 40- 85% Cu rich Co-Cu solid solutions ||l| and 
variation of magnetic moments for Co rich alloys at room 
temperature Recently, the focus of attention has 

been on giant magneto-resistance (GMR) studies. Sev- 
eral experimental groups have been carrying studies of 
GMR in this system. Most remarkable works include in- 
vestigation of GMR in bulk CoajCui-a; alloys for x = 
0.05-0.2 [^,||, in heterogeneous thin film CoCu alloys 
[^,0 and in as-grown, epitaxial Co-Cu alloy layers 
Other magnetic studies on this system include magnetic 
anisotropy of Co on Cu(llO) ^ and of Co-Cu ultra thin 
films [Q. An altogether different type of work by Chil- 
dress and Chien pd| ] is, to our knowledge, the only one of 
its kind which involved low-temperature magnetic stud- 
ies of Co^^Cui-a; alloys. This work is worth mentioning 
because it revealed very interesting features of the mag- 
netic phases over the entire concentration range from x 
= to x= 0.8. It indicated a low-temperature spin glass 
phase upto x = 0.23 ; a low-temperature re-entrant spin- 
glass or mixed phase, stable in the range 0.24 < x < 0.40 
and the normal low-temperature long-ranged random fer- 
romagnetic phase stable beyond x = 0.4. Mandal and 
Ghatak have reported calculations on the periodic 
Anderson model for a binary alloy, working within the 
Hartree-Fock, coupled with the Virtual Crystal Approx- 
imation (VGA). They have been able to produce some 
of the features observed experimentally. The model used 
by these authors have many empirical fitted parameters. 
The existence of a spin-glass and, at higher concentra- 
tions, a re-entrant mixed phase has been known from sev- 
eral mean- field approaches |13 3]. To our knowledge, the 



system has never been investigated from a first-principles 
electronic structure point of view. The aim of this com- 



munication is to present such a calculation of magnetic 
properties of the fee, ferromagnetic phase of Co^Cui-a; 
alloys in the concentration range from 40% to 80% of 
Co. We have used the fully self-consistent, tight-binding 
linearized muffin-tin orbitals augmented space recursion 
(TB-LMTO-ASR) technique as the basis of our calcula- 
tions. We have restricted ourselves to the ferromagnetic 
phase only leaving the spin-glass and mixed phases for a 
later study. Our studies involved calculation of magnetic 
moments. Curie temperatures and the spin susceptibili- 
ties. 



II. THEORETICAL DETAILS 
A. The magnetic phases 

Description of magnetic phases within the local spin 
density approximation (LSDA) involves the study of the 
evolution of local magnetic moments in the vicinity of ion 
cores because of the distribution of the valence electron 
charge. Each lattice site in the face centered cubic struc- 
ture is occupied by an ion core : in our case randomly 
by either Co or Gu. We shall associate a cell or a sphere 
with each ion core and assume that the charge contained 
in the sphere belongs to that ion core alone. Ideally 
such cells or spheres should not overlap. In the tradi- 
tional Kohn-Korringa-Rostocker (KKR) method this is 
certainly so. However, in the atomic sphere approxima- 
tion (ASA) which we shall use in our TB-LMTO ver- 
sion, this division of space is to a certain extent arbitrary. 
Within these cells the valence electrons carrying spin a 
sees a binary random spin-dependent potential V^(r), 
where A = Go or Cu and a = f or j. 

The charge density within the cells can be obtained 
from the partially averaged Green functions : 
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~p,{r) = -(1/^)3™^ / [x « Gfr E) » 



(1) 



where < G^°'"(r, r, £;) > and < G^l'" [r,r, E) > 
are partially averaged Green functions with the site r 
occupied by a Co or Cu ion core potential corresponding 
to spin (T. 

For the random ferromagnetic phase we proceed as fol- 
lows : we consider all cells to be identical in that they all 
carry identical average charge densities. We shall borrow 
the notation of Andersen etal [|l5| to write functions like 
f{fR) which are equal to /(r) when r lies in the atomic 
sphere labelled by R and is zero outside. The ferromag- 
netic charge densities are defined as : 



R 

P2{r) = ^PiirR) 



R 



The magnetic moment per cell (atom) is then defined 
by: 



m^{l/N) J d3?:[pi(r) - P2W] 

= {l/N)J2 I d'^rlp^im) - PiiTR) 
= (l/iV)V/ mfl(rfl) 

„ Jr<S 



Since all cells are identical, the above calculation need 
be done only in one typical cell. Within the TB-LMTO- 
ASA the cells are replaced by inflated atomic spheres 
and the remaining interstitial is neglected. The problem 
is then one of a binary alloy with an almost non-magnetic 
charge density due to the Cu ion cores and a magnetic 
one due to the Co ones. The averaging is done over con- 
figurations of the random alloy. 



B. The configuration averaging and TBLMTO-ASR 

For random alloys we extract physical properties using 
configuration averaging which means that the physical 
quantities involved are the configuration averaged quan- 
tities. A powerful technique of carryingout this averag- 
ing is the augmented space recursion [|l6| . The method 
allows us to go well beyond the traditional single site 
coherent potential approximations and has been applied 
successfully to a wide variety of systems 
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The 

convergence of the ASR has been established recently 
pl[ , so that any approximation we impose on the recur- 
sion is controlled by tolerance limits preset by us. ASR 
coupled with TB-LMTO-ASA has been proved to be a 
very powerful technique in predicting material properties. 



The TB-LMTO-ASR has been described in great de- 
tail earlier ||l6|- |l9| , ^ . We shall refer the reader to the 
referenced monograph for technical details. Here we quote 
only the main results. 

Our starting point is the TB-LMTO Hamiltonian in 
atomic sphere approximation (ASA) : 



E„ 



hi 



where, h = C 



E, 



AI/25AI/2. 



C , and A are diagonal matrices in angular momentum 
space and are th e potential parameters of the TB-LMTO 
technique and S is the structure matrix which is sparse 
in the most tight-binding representation. 

Let us now look at the transformation of the Hamil- 
tonian in the full augmented space when the potential 
parameters have homogeneous binary random distribu- 
tions. Such random variation may be described by ran- 
dom site occupation variables tie. which take values 1 or 
according to whether the site R is occupied by an A-type 
or B-type of atom and have probability densities 

p[nR) = X 6{nR - 1) + (1 - x) S{nR) 

The details of construction of the augmented space and 
the effective Hamiltonian on it has been described in de- 
tail in our earlier communications []l8| , ^3|j2^ ] . Wc shall 
quote here only the final result and refer the readers to 
the referenced works . The trick involves first in replacing 
each random parameter P in the Hamiltonian by expres- 
sions of the type : 

Paur + PB{l-nR) = Pb + [PA-PB)nR (2) 

and then replacing each factor ur by an operator AAr 
which acts on the configuration space of n^. For a binary 
distribution of hr the rank of this space is two. We shall 
construct a basis in this space which we shall designate 

{I T) 1 I i)}- For the binary distribution this operator 



Mr = xvf + {i~x)vf + ^x{i - x) {r« + r^f} 

and the above equation becomes : 

= A{P)I + B{P)Vf + F(P){r^f + T^f} 



P 



where 



A(P) 
B(P) 
F(P) 



xPa + {1-x)Pb 
(1 - 2x) {Pa - Pb) 



^x(l-x) (Pa-Pb) 
We now proceed as follows : 

{E - H)-^ ^{E - C - A^/^^A^/^)-! 

'E-C 



A-1/2 



S 



A-1/2 
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All factors except S have binary random distributions. 
Using the procedure described above we may now con- 
vert the above equation into augmented space. The aug- 
mented space theorem gives P3| 

« Gi^L,flL(£^) > = {R,L,m\(^EI - \R,L,{iD}) 
First note that : 

A-'/^\R, L, {$}) - A(A-i/2)|i?, L, {$}) 

+ F(AV2)|i?,L,{i?}) = 11} (4) 

1 /2 

and if we define [A(l/A)] ' |1} as then this latter 
ket is normalized. A little algebra then gives : 

^Grl,rl{.E)-> ^ {l\[E-A + B + F-S]-^\l) (5) 

where 

A = {A(C/A)/A(1/A)}I®Z®X 

RL 
RL 

S - EE{A(l/A)-^/^}W'L'{A(l/A)V2}r^^, 

RL R'L' 

Tll"»X (6) 

and J = Ja + Jb + Jf where the operators on 
the right hand side of these equations are given by 

Ja = {A{{C - E^)/A) /A{l/A)} I(E)I(E)I 

Jb - {B((C-£;,)/A)/A(l/A)} E^i^^-^i^^f 

RL 

Jf = {F{{C^E,)/A)/A{1/A)}J2'Pr'»Tl 

RL 

This form is suitable for recursion. We start the recur- 
sion with the state |1) which is a mixture of the configu- 
rations {0} and {R} and carry out recursion as usual. 

The initial TB-LMTO potential parameters are ob- 
tained from suitable guess potentials as described in the 
article by Andersen etal [|l5|. In subsequent iterations 
the potentials parameters are obtained from the solution 
of the Kohn-Sham equation 

{-^\7^ + V''^ -E]rArR,E) - (7) 
where, 

V^irR) = VZeirR) + V;::ArR) + V,7{rR) + V,nad 

(8) 



here v refers to the species of atom sitting at R and 
a the spin component. The electronic position within 
the atomic sphere centered at R is given by tr =r — R. 
The core potentials are obtained from atomic calculations 
and are available for most atoms. The treatment of the 
Madelung potential in a random alloy has always raised 
problems. We adopt a procedure suggested by Drchal 
etal |2j] and regularly used in CPA calculations within 
the TB-LMTO. We choose the atomic sphere radii of 
the components in such a way that they preserve the 
total volume on the average and the individual atomic 
spheres are almost neutral. This ensures that total charge 
is conserved, but each atomic sphere carries no excess 
charge. However, we are careful that such a choice does 
not violate the overlap criterion of Andersen and Jepsen 
JlSf . In the ASR-LSDA self- consistency loop, charge 
transfer takes place between these spheres ; however, at 
the end of the self-consistency iterations, the spheres are 
approximately neutral and hence do not contribute to a 
Madelung energy. This prescription is to an extent ad 
hoc, and there is no guarantee in general that we will 
be able to find such atomic sphere radii. However, the 
procedure has proved rather successful in many earlier 
CPA |2^] and ASR calculations on magnetic alloys 
and we shall adopt it here. As in the CPA calculations we 
iterate until the total energy and moments of the charge 
density converge. In this sense our calculations are self- 
consistent in the LSDA sense. 



C. Curie temperature 

We have estimated the Curie temperature by two 
different methods-Bragg Williams approximation and 
Mohn-Wolfharth model ||2^. Both the procedures are 
described below: 

1. Bragg- Williams approximation: 

The estimates of the Curie temperature were ob- 
tained from the magnetic pair energies ||2^ The 
pair energies are defined as follows : At two sites 
labeled r and r' in a completely random paramag- 
netic background, we replace the potential by that 
of either the up-spin ferromagnetic Co or the down- 
spin one. The Green function of this system we 
shall denote by : G^^'*^*^ (r, r, E), a being the spin 
type at the site r (either f or j) and a' that at the 
site r' . The pair energy is defined as 

E{R) = dE E[~{l/n)^m{Gl'i:^\r,r,E) 

J —OO 

+ G'll^'Hr,r,E)-G^^f'{r,r,E) 

~ Gl''^^\T,r,E))\ (9) 

Here R — r — r' . We may either estimate the 
above directly, or to be more accurate we may use 



3 



the orbital peeling method of Burke |gj . The lat- 
ter is an extension of the recursion method, where 
small differences of large energies (as in the def- 
inition of the pair energy) are obtained directly 
and accurately from the recursion continued frac- 
tion coefficients. Note that we have assumed that 
the dominant contribution to the pair energy comes 
from the band contribution and the rest approxi- 
mately cancel out. The simplest Bragg- Williams 
estimate of the Curie temperature is 



Tr 



{1~x)E{Q)/kb 



where 

£;(0) = E{q = 0) and E{q) = Y,Rexp{iq.R)E{R). 
Since the pair energy is short-ranged, a reasonable 
estimate of E{0) is J2n<3 ^i^n) where Rn is 
the n*''-nearest 

neighbour vectors and Z„ is the number of n"*- 
nearest neighbours. The Bragg- Williams approach 
overestimates the Curie temperature and its gener- 
alization, the cluster variation method, yields bet- 
ter quantitative estimates. We have restricted our- 
selves to the Bragg- Williams nearest neighbour pair 
energy approximation. 

2. Mohn-Wohlfarth model(MW): 

In this model [p^ the Curie temperature is given 
by: 



rp2 



Tsi 



where, 

T^ is the stoner Curie temperature given by, 

nEF)r^N{E)[IL)dE = 1 

I(Ej') is the Stoner parameter obtained from the 
earlier calculations , N(E) is the density of 
states per atom per spin |30(| and f is the Fermi 
distribution function. 

TsF is the spin fluctuation temperature given by. 



Xo is the exchange enhanced spin susceptibility at 
equilibrium and m is the magnetic moment per 
atom. 

Xo is calculated using the relation by Mohn and 
Gersdorf 



I is the stoner parameter and 'N^{Ef) and '^'^{Ef) 
are the spin-up and spin-down density of states per 
atom. 



III. COMPUTATIONAL DETAILS 

For the calculation of the component projected av- 
eraged density of states of the ferromagnetic phase we 
have used a real space cluster of 400 atoms and an aug- 
mented space shell upto the sixth nearest neighbour from 
the starting state. Eight pairs of recursion coefficients 
were determined exactly and the continued fraction ter- 
minated by the analytic terminator due to Luchini and 
Nex In a paper Ghosh etal have shown the 

convergence of the related integrated quantities, like the 
Fermi energy, the band energy, the magnetic moments 
and the charge densities, within the augmented space re- 
cursion. The convergence tests suggested by the authors 
were carried out to prescribed accuracies. We noted that 
at least eight pairs of recursion coefRcients were neces- 
sary to provide Fermi energies and magnetic moments 
to required accuracies. We have reduced the computa- 
tional burden of the recursion in the full augmented space 
by using the local symmetries of the augmented space 
to reduce the effective rank of the invariant subspace in 
which the recursion is confined JT6t and using the seed re- 
cursion methodology ||2^ with fifteen energy seed points 
uniformly across the spectrum. 

We have chosen the Wigner-Seitz radii of the two con- 
stituent atoms Co and Cu in such a way that the aver- 
age volume occupied by the atoms is conserved. Within 
this constraint we have varied the radii so that the final 
configuration has neutral spheres. This eliminates the 
necessity to include the averaged Madelung Energy part 
in the total energy of the alloy. The definition and com- 
putation of the Madelung Energy in a random alloy had 
faced controversy in recent literature and to this date no 
satisfactory resolution of the problem exists. Simultane- 
ously we have made sure that the sphere overlap remains 
within the 15% limit prescribed by Andersen. 

The calculations have been made self-consistent in the 
LSDA sense, that is, at each stage the averaged charge 
densities are calculated from the augmented space re- 
cursion and the new potential is generated by the usual 
LSDA techniques. This self-consistency cycle was con- 
verged in both total energy and charge to errors of the or- 
der 10~^. We have also minimized the total energy with 
respect to the lattice constant. The quoted results are 
those for the minimum configuration. No short ranged 
order due to chemical clustering has been taken into ac- 
count in these calculations, nor any lattice distortions 
due to the size differences between the two constituents. 



Xo = 



2^l% \2m{EF) ' 2NI{Ef) 



- I 
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IV. RESULTS AND DISCUSSION 

Figure 1(a) and 1(b) show the partial density of states 
at Co and at Cu sites respectively for various concentra- 
tions of Cu. It is evident from the figures that there is 
hardly any difference in the relative heights and shapes 
in the majority and minority spin bands at Cu site for 
various concentrations. This is reflected in the local mag- 
netic moment shown in Figure 2 where Cu sites have a 
very small amount of induced moment. For the minor- 
ity spin partial densities on Co, the structure below the 
Fermi level do not change much with the Co concentra- 
tion. The peak at around 0.0 Ryd. grows with decreas- 
ing Co concentration. However, this structure is above 
Fermi level and does not contribute to the magnetic mo- 
ment. For the majority spin partial densities on Co, the 
peak around -0.4 Ryd. doesn't change with concentration 
while the one around -0.3 Ryd. grows with concentration. 
The most remarkable change occurs in the peak around 
-0.2 Ryd. whose height reduces with the concentration 
and finally the peak around -0.3 Ryd. contributes more 
with increasing concentration. As a result, the partial 
magnetic moment on Co increases only slowly with con- 
centration and the variation of average magnetic moment 
is linear as shown in Figure 2. 

Figure 3 shows the variation of average magnetic mo- 
ment in emu/gm unit and the agreement with experimen- 
tal results |1^ is pretty good qualitatively. Both CPA 
and ASR moments agree well too. But, the agreement 
with the results from model calculation JT2| is far from 
being reasonable. Their calculations predict a monotonic 
fall of average magnetic moment and finally vanishing 
of magnetic moment around 40% of Co concentration. 
This has neither been observed in experiments pl nor 
in our calculations with two widely used ab-initio meth- 
ods. This discrepancy is due to the choice of technique 
in case of model calculations. They used Hartree-Fock 
and VCA for their calculations. Both the approximations 
have their own limitations in predicting properties of real 
materials. The very important electron-electron correla- 
tion effect is absent within the framework of Hartree- 
Fock approximation while in LSDA this effect has been 
incorporated. The limitation of VCA is that it has been 
proved to be successful only in cases of weak scattering 
alloys where each particle sees nearly the same average 
perturbation field. ASR on the other hand is a very pow- 
erful and accurate technique in the sense that it is an ex- 
act method involving no mean-field like approximation 
for configuration averaging. This is the reason for bet- 
ter agreement of our calculations with the experimental 
results. 

Figure 4 shows the partial and average spin susceptibil- 
ities of the system. The partial susceptibility of Cu varies 
smoothly and finally becomes almost saturated while we 
find little oscillations in case of Co particularly in the 
range of 60%-80% of Co. This nature is reflected in the 
average susceptibility. This is basically due to oscilla- 



tions in the up density of states at Fermi energy at Co 
sites(See Table I). However, the oscillation in the Co as 
well as average susceptibilities are very small, eventually 
leading to saturation. 

Figure 5 shows the variation of Curie temperature with 
Co concentration using Bragg Williams (BW) approxi- 
mation and MW model. In BW approximation, the vari- 
ation of Curie temperature is smooth though the values 
are overestimated. This is expected as BW approxima- 
tion has a tendency to overestimate j3^. This overesti- 
mation is reduced using MW model but the qualitative 
agreement between the two methods is not obtained. In 
the concentration regime 40%-60% of Co, Curie temper- 
ature oscillates in MW model and beyond it qualitative 
agreement with BW calculations is achieved. This fluc- 
tuating behaviour obtained in MW calculations is due to 
the oscillations in Spin fluctuation temperature(See Ta- 
ble II) because Stoner Curie temperature has a system- 
atic variation with Co concentration. This oscillation in 
spin fluctuation temperature can be explained as follows: 
for the concentration 40%-60% of Co the partial spin sus- 
ceptibilities of both Co and Cu undergo variations which 
is reflected in the spin fluctuation temperature while for 
60%-8Q% of Co the susceptibilities of both the compo- 
nents remain almost constant except for a small variation 
at around 80% of Co. That's why the spin fluctuation 
temperature and hence the Curie temperature follow the 
same trend as obtained in BW calculations for Co concen- 
tration beyond 60%. The Stoner Curie temperature(See 
Table II)is very high as expected because LSDA overesti- 
mates the binding. We have also tried to characterize the 
system for the concentration regime of our investigation 
following the prescription by Mohn and Wolfharth |2^] 
which characterizes a system from the values of Curie 
temperature and spin fluctuation temperature. Accord- 
ing to the prescription the quantity tc — determines 

c 

the agency playing the dominant role in finding out Curie 
temperature and hence the effect of spin fluctuations in 
determining the magnetic properties of the system. The 
systems with tc <0.5 are called fluctuation systems in 
which description of magnetism should involve spin fluc- 
tuation effects and systems with tc >0.5 are called Stoner 
systems in which single particle excitations play the dom- 
inant role. In our case, we see (See Table II) that for the 
entire range of concentrations tc remains less than 0.5 
which means that the description of magnetic properties 
should have the effect of spin ffuctuations. 

V. CONCLUSIONS 

We have calculated for the first time, the electronic 
structure and magnetic properties of Co^Cui^a; alloys 
using fully self-consistent first principles electronic struc- 
ture technique and obtained results which agree reason- 
ably well with the experimental results. This clearly 
shows that ASR coupled with TBLMTO is a powerful 
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technique in describing electronic structure and magnetic 
properties of binary alloys. Our results clearly show the 
flaws in the theoretical results based upon model calcula- 
tions. Finally, this study clearly presents the fact that in 
CoCu alloy spin fluctuation effect is very much present in 
the ferromagnetic phase and a more accurate description 
of the magnetic properties of the system should involve 
this eff'ect. 
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FIGURE CAPTIONS 



VI. TABLES 



Figure 1(a) shows the partial densities of states of 

Co at various Co concentrations. The concentrations are 
indicated in the inset. Dashed curves show the results 
for the minority and full curve the majority spin states. 
Fermi energies are indicated by vertical dashed lines with 

each figure. 

Figure 1(b) shows the partial densities of states of 
Cu at various Co concentrations. The concentrations are 
indicated in the inset. Dashed curves show the results 

for the minority and full curve the majority spin states. 
Fermi energies are indicated by vertical dashed lines with 
each figure. 

Figure 2 shows The partial magnetic moments on Co 
and Cu sites and the average magnetic moment as a func- 
tion of Co concentration. The solid and dotted lines refer 
to the augmented space calculations while the squares, 
the crosses and diamonds to the CPA. All these values 
are in units of Bohr-magneton. 

Figure 3 shows the average magnetization in emu/gm 
as a function of Co concentration at OK obtained from 
the theoretical estimates. The solid line refers to the ASR 
calculations while the diamonds to CPA. 

Figure 4 shows the average and partial spin 
susceptibilities (2^) in Ryd~^atom~^ with Co concen- 
tration. 

Figure 5 shows the Curie temperature as a function 

of Co concentration. The solid line corresponds to the 
results using BW approximation while the squares corre- 
spond to the results with MW model. 
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8.17 



TABLE I. shows the density of states (both up and down) 
at fermi level for both the components and the partial and 
average susceptibilities with Co concentration 





ryBW 

c 


TsF 


^ C 


rpMW 


tc 


0.4 


1042 


1169 


7007 


1138 


0.16 


0.45 


1150 


1089 


7095 


1064 


0.15 


0.5 


1394 


1134 


7186 


1107 


0.16 


0.6 


1667 


1016 


7356 


998 


0.14 


0.7 


1911 


1260 


7530 


1227 


0.17 


0.8 


1983 


1728 


7709 


l()i8 


0.22 




TABLE 




II. 






shows 



the values of BW Curie temperature(Tf MW Curie 
temperature(T^^),spin fluctuation temperatures(TsF) and 
Stoner Curie temperature(Tf ) with Co concentration 
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Figure 1(a) 
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Figure 1 (b) 
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Figure 5 




